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RELATIVISTIC CHARGED FLUID FLOW 


II: Generalized Larmor and Helmholtz Theorems 

by 

Lawrence A. Schmid 
ABSTRACT 

The covariant Larmor theorem is derived within the framework of Special 
Relativity for an inviscid charged fluid in the presence of an arbitrary electro- 
magnetic field and a gravitational field which is described by a simple scalar 
potential. The flow need not be adiabatic and the fluid need not be barotropic. 
The covariant form of the Larmor theorem, which is an antisymmetric tensor 
equation, is equivalent to two 3-vector equations, one of which is the fluid equa- 
tion of motion. The other is just the familiar 3-vector statement of the Larmor 
theorem except for an extra precession that is thermal in origin. The relativis- 
tic Helmholtz equation is derived from the covariant Larmor theorem. The vor- 
ticity is defined as the curl of the usual canonical momentum of a charged par- 
ticle in an electromagnetic field except that the particle rest mass is a variable 
that includes contributions that are proportional to the gravitational potential and 
the specific enthalpy of the fluid. For this definition of the vorticity it is found 
that the vorticity flux diffuses through the fluid except when the flow is isentropic. 
It is shown, however, that when an appropriate thermal 4-potential is included in 
the definition of the canonical momentum, the vorticity flux remains frozen in 
the fluid even when viscosity is present and the flow is nonadiabatic. A micro- 
scopic interpretation of this generalized Helmholtz equation is given. 
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RELATIVISTIC CHARGED FLUID FLOW 


It: Generalized Larmor and Helmholtz Theorems 
I. INTRODUCTION 

In the preceding paper,* henceforth referred to as I, it was shown that the 
entropy-dependent force terms in the relativistic equation of motion could be 
expressed in terms of an antisymmetric entropy force tensor whose dynamical 
role is analogous to that of the electromagnetic field tensor. This analogy is 
most clearly evident in the relativistic statement of the Larmor theorem for a 
charged inviscid fluid, which is derived in Section II of the present paper. 

Use of the entropy force tensor also makes possible the derivation of the 
relativistic generalization of the Helmholtz equation for the time dependence of 
the fluid vortieity. This is carried out in Section HI. In the form of the Helm- 
holtz equation derived in Section III, the vortieity is defined as the curl of the 
canonical particle 4-momentum /xv 1 + (q/c)A j where /uis the particle rest mass, 
including contributions from the gravitational and thermal energy, q is the parti- 
cle charge, and A j is the electromagnetic 4-vector potential in Gaussian units. 

As in I, v> denotes the particle 4-velocity and c, of course, is the speed of light. 

The form of the Helmholtz equation in terms of a vortieity defined as the 
curl of /j-v 1 + (q/c) A j has the disadvantage that the vortieity so defined is in 
general not simply carried along with the fluid flow, but rather diffuses across 
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flow lines. In Section IV, however, it is shown that it is possible to define the 
vorticity in such a way that no diffusion occurs. The vorticity having this prop- 
erty is defined as the curl of a generalized canonical momentum that includes 
not only the contribution (q/c) A' from the electromagnetic field, but also an 
analogous contribution that accounts for viscosity as well as for the entropy- 
dependent forces. 

Notation 

The notation employed will be the same as that used in I with the exception 
of a modification in the notation for temperature. Inasmuch as we have com- 
mitted ourselves to a 4-vector description of temperature, it would be desirable 
to introduce a notation that is analogous to that for particle mass and density. 

In the fluid rest frame these are designated by m and p respectively, and in the 

observer's frame by m* and p*. Thus temperature in the fluid rest frame will 

° * 
be designated by T (instead of T as in I), and in the observer’s frame by T (in- 
stead of T). Equation (3.12) of I will now be written 

V = T vi/c = T*(l, v/c) (1.1a) 

where 

T* = TT, (1.1b) 

and 

r = (1 - v 2 /c 2 y 1/2 . (1.1c) 

Similarly, the heat reservoir will be described by the 4-vector P where 

T i = T k ^/c = T* R d, .„/c), (1.2a) 
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(1.2b) 


t = r t . 

R R R 


T r = (1 - v2/ C 2)-l/2, 


where obviously v R is the 3-velocity of the reservoir, T R is its temperature in 
its own rest-frame, and T* is its temperature in the observer's frame. The 
4-vector must satisfy the reversibility condition given in (4.27a) of I: 

R 

V. T j = V. Ti = cT. (1.3) 

i R J 

As in I, the antisymmetric electromagnetic field tensor F jk is expressed in 
terms of the electric and magnetic field intensities E and B as follows: 


(F 10 , F 2 °, F 30 ) = (E x> E y , E z ) = E; 


(F 23 , F 31 , F 12 ) 3 _ (B X , B y , B z ) = - B 


If A j is the electromagnetic 4-potential, then 


F jk = 3i A k - d k A j . 


In terms of the scalar potential A° and the 3-vector potential A 


AJ = (A 0 , A). 


From (1.4) - (1.6) it follows that 


E = - VA° - 'dA/cBt ; 


(1.7a) 


B = V x A . 


(1.7b) 
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We shall have occasion to decompose the tensor F jk into the two antisym- 
metric tensors F jk and F jk which are its parallel and normal projections with 
respect to the 4-velocity v j . The tensor F jk is characterized by the condition 

F jk v k = F jk v k . (1.8a) 

Thus the electromagnetic ’force calculated using F jk is the same as that calcu- 
lated from F jk . The condition (1.8a) alone does not suffice to specify F ik uniquely. 
It must be augmented by the condition that in the fluid rest frame the space-space 
components of F jk vanish, i.e. 

(F 23 , F 31 , F 1 2 ) v _ Q = 0. (1.8b) 

From (1.8a) it follows that 

(F 10 , F 20 , F 30 ) v = 0 = E, (1.8c) 

O 

where E is the electric field intensity in the fluid rest frame. It is easily veri- 
fied that F J k must have the following form: 


or 

F jk = (F J ^v^v k + v j v^F^ k )/c 2 ; 

(1.9a) 


(F 10 , F 20 , F 30 ) = T E , 

(1.9b) 

where 

(F 23 , F 31 , F 12 ) = T E x v/c, 

(1.9c) 


E = T [(E-Bx v/c) - (E • v ) v/c 2 ] . 

(1.9d) 
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Note that E is just the effective electric field intensity that appears in the rela- 
tivistic form of Ohm's Law. 1 

The antisymmetric tensor F i k is given by 

F' k = F jk - F jk . (1.10a) 


From (1.4) and (1.9) we find 


(F 10 , F 20 , F 30 ) = r B X v/c, 

(1.10b) 

(F 23 , F 31 , F 12 ) = - T B, 

(1.10c) 

where 


B = r[(B + Ex v/c) - ( B • v ) v/c 2 ] . 

(l.lOd) 

In the fluid rest frame we have 


(F 10 , F 20 , F 30 ) v = 0 = 0, 

(1.11a) 

(F 23 , F 31 , F 12 ) v = 0 = -B. 

(1.11b) 


Thus F jk is the tensor to be associated with effects such as Larmor precession 

O O 

that in the fluid rest frame depend only on B and not on E. We shall refer to E 

and B as the effective electric and magnetic field intensities. 

% 

Fluid Rotation Tensor 

As a preliminary to deriving the relativistic Larmor theorem in the next 
section, we shall now discuss the rotation tensor for the fluid. First we note that 
if X> = (0, r ) is the space-time displacement between two points that are observed 
simultaneously (X° = cAt = 0), and if one point is rotating about the other (the 
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origin) with angular velocity Q as seen by an observer for whom the origin has 
the velocity v , then 

dr/dt = £2 xr. (1.12) 

The left side of (1.12) may be put into covariant form by noting that dr = dt/P 
where dT is the proper time interval in the rest frame of the origin. Thus 
Tdr/dt = dr /dr. If we now introduce an antisymmetric tensor whose space- 
space components are given by 


(n 23 , n 31 , n 12 > = ra, (i.i3) 

then (1.12) can be written in the following covariant form: 

dXVdr=ni k X k . (1.14) 

The space part of (1.14) is just (1.12). It follows from the antisymmetry of H )k 
that 

0 = X. dXJ/d r = IdCXj Xj)/dr = idr 2 /dr. (1.15) 

This is the necessary condition that the length of the radius vector of a rotating 
point be constant. 

It is well known that the non- relativistic definition of Q in the case of a fluid 
is given by fi = ( 1/2) Vx v . The obvious covariant generalization of this is 

Q jk = - I (3i v k _ 3 ^ v j ); (1.16a) 
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or 


(ft 10 , n 20 , fi 30 ) = [3<rv>/3t + V(r c 2 )]/2c, (l.ieb) 

(f) 23 , n 31 , n 12 ) = i v x (Tv) = ro. (i.i6c) 

The last equality in (1.16c) constitutes the definition of the angular velocity Q 
in the relativistic case. The factor T has been inserted into this definition in 
order to give (1.16c) the same form as (1.13). 

In the same way that we decomposed F jk , we shall decompose fi jk into two 
antisymmetric tensors Q ik and Q ik which are the parallel and normal projec- 
tions with respect to v j . The tensor fl jk is defined by the conditions 

fi jk v k = n jk V k = J d V j /dr (1.17a) 

and 

(fi 23 , fi 31 , 0 12 ) v=o = 0. (1.17b) 


The last step of (1.17a) follows from (1.16a) above and (1.5), (1.6a) of I. It is 
easily verified that Q jk has the following form: 

Qjk _ (fii^ V/ £ v k + v 1 v^fA)/c 2 

= [(d vVdr) v k - (d v k /dr) v 1 ] /2c 2 ; (1.18a) 

or 

(fi 10 , f) 20 , Q 30 ) = (r 2 /2c) d y/dr, (1.18b) 

(n 23 , n 31 , q 12 ) = r n T (i.i8c) 
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where 


Q t = - (r 2 /2c 2 ) v x Dv (1.18d) 

where D = P ~ 1 (cl/dr) is, as defined in (1.6b) of I, the substantial time differentia- 
tion operator in the observer’s frame. Note that ffi T is just the well-known 
Thomas precession. 2 

The antisymmetric tensor Q > k is given by 

6) k = Q i k - f}i k . (1.19) 

From (1.17a) and (1.19) 

5i k v k = 0. (1.20) 

Any antisymmetric tensor satisfying the condition (1.20) has only three degrees 
of freedom and can be written in the form 

(fi 10 , Q 20 , Q 30 ) = (T/c) V X Q, (1.21a) 

(h 23 , q 31 , n 12 ) = r q, (i.2ib) 

where (1.21b) constitutes the definition of the 3-vector Q in terms of the space- 
space components of f2i k . From (1.19), (1.16c), (1.18c), and (1.21b) it follows 
that 

n = q + n T . (f- 22 ) 

In the next section we shall see that Q is the fluid rotation that is produced by 
the magnetic and thermal fields plus the residual rotation that has been retained 
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from the initial conditions of the fluid. Equation (1.22) simply says that in order 
to arrive at the total rotation £2, we must also include the Thomas precession £2 T . 

Entropy Force Tensor 

It will be necessary to carry out the same kind of decomposition for the 
entropy force tensor 0’ k which in (5.12) of I was defined as 

0 jk = T^ 3 k s - T k 3 j s. (1.23) 

From (4.26a) of I we have, making the notation change noted in writing (1.1), 

t r = ( T /c) ( vi + w, )> (1.24a) 

or 

= T* [(1 +v • w/c 2 ), (v +w)/c] (1.24b) 

where T* = TT is the temperature in the observer’s frame and v + w = v R is the 
reservoir 3-velocity. The 3-velocity w is obviously the drift velocity with re- 
spect to the fluid. Using (1.24) and (1.23), we arrive at the following relations: 


0 jk = (T/c) [(v j +wJ) 3 k s - (v k + w k ) 31 s] ; (1.25a) 

(0 10 , 0 20 , 0 30 ) = T* Vs + (TVc 2 ) [(v + w) 3s/3t + (v • w)Vs] ; (1.25b) 

(0 23 , 0 31 , 0 12 ) = (T* V s) x (v+w)/c. (1.25c) 

The tensor © jk is defined by the conditions 

0 jk v k = ©■> k v k = s - cTBi s (1.26) 
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and 


(0 23 , 0 31 , 0 12 ) v=o = 0, (1-27) 

where the second step of (1.26) follows from (5.13) of I. Thus 

@jk _ + v j v^0^ k ),/c 2 (1.28a) 

= (T/c) [ (wd s/c 2 - s) v k - (w k s/c 2 - 3 k s) ] ; (1.28b) 

(0 10 , 0 20 , 0 30 ) = r 0, (1.28c) 

(0 23 , 0 31 , 0 12 ) = r© X v/c, (1.28d) 

where 

0 = TVs + (me 2 ) -1 [(mT*s)w+ (mT3s/<3t-mT*sv •w/c 2 )v], (1.28e) 

The antisymmetric tensor 0 jk is given by 

©jk = 0jk _ 0ik . (1.29) 

From (1.26) and (1.29) it follows that 

0 jk v k = 0. (1.30) 

Using (1.25) and (1.28) in (1.29) we find 

0j k - (T/c) I w J ( B k s - v k s/c 2 ) - w k ( cd s - v k s/c 2 )] ; (1.31a) 

(0 10 , 0 20 , 0 30 ) - - TO x v/c, (1.31b) 

(0 23 , S 31 , 0 12 ) - TO, (1.31c) 
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where 


0 = (T V s + T* s v/c 2 ) x w/c . (1.31d) 

It is evident that 0 jk can be written in terms of the 3-veetor 0 and 0 as follows: 

(0 10 , 0 20 , 0 30 ) = T (0 - 0 x v/c), (1.32a) 

(0 23 , 0 31 , 0 12 ) = T (0 + 0 X v/c). (1.32b) 

Note that in the above relations the relative velocity of the reservoir w j or 
w may be replaced by the thermal force <p j or 9 by means of the relation (4.23) 
of I: 

cp j = (mTs/c 2 ) w' ; (1.33a) 

9 = (mTs/c 2 ) w. (1.33b) 


Gravitation 

It would be desirable to include the effects of gravitation to the extent that 
they are of importance in practical magnetofluid-dynamical problems. For this 
purpose General Relativity is not needed. A simple flat-space covariant scalar 
theory will suffice. A variety of such theories have been developed . 3 The first 
and simplest of these, and the one that we shall employ, is Nordstrom’s theory . 4 
This theory requires that the rest mass /lx of a particle be exponentially dependent 
on the scalar gravitational potential G, i.e /x = m exp (G/c 2 ) where m is a constant. 
The necessity for this dependence follows from the equation of motion (in the 
absence of all fields but the gravitational one): d^v^clr = /x^G. Contraction 
with v j gives d/u/dr = (/x/c 2 ) d G/dr, and integrating this yields /lx = m exp (G/c 2 ) 
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where m is the integration constant. Note that if the equation of motion had been 
taken as d(/xv j )/dr = m3 j G, where the coupling is given by the constant mass m 
rather than the variable mass /x, then we would have been led to the linear rela- 
tion /-x = m (1 + G/c 2 ) rather than the exponential one given above. The difference 
between these two relations is of order m (G/c 2 ) 2 , which for our purposes is com- 
pletely negligible. Thus we shall adopt the simpler linear relation, although in 
principle the exponential one is preferable because it is consistent with the Equi- 
valence Principle in that the same mass /x that describes the inertial properties 
of the particle also describes the strength of its coupling to the gravitational 
field. 

Referring to (5.11) of I, we recall that the variable particle mass m = 
m (1 + h/c 2 ) both determined the inertial properties of the particle and also 
served as a thermal potential function. We may now define a total particle rest 
mass /x that includes both the gravitational and thermal contributions as follows: 

/x = m(l + G/c 2 + h/c 2 ), (1-34) 

with the result that 

cM /xc 2 - mcdG + mB-'h. (1.35) 

Thus the 4-gradient of the variable particle rest energy /x c 2 gives both the gravi- 
tational force and the enthalpy-dependent thermal force. The fluid equation of 
motion including gravitation is just (5.17) of I with m replaced by /x: 

d(/xv j )/d r = cb (fic 2 ) + (m/c) 0 ik v k 

- (\ s > k )/p + (q/c) F jk v k . (1.36) 
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Because no assumption has been made regarding the equation of state of the 
fluid, we are not restricting ourselves to the most commonly considered special 
case of a barotropic fluid. 

In the next section we shall find that the forces arising from G and h are 
represented in the relativistic Larmor theorem by means of the tensor G J k which 
is defined as follows: 


G jk = (c/m) [(Bi /x) v k - ( B k /x) v j ] 


= (B - 1 G + B* h) v k /c - ( B k G + B k h)vVc; (1.37a) 

(G 10 , G 20 , G 30 ) = - r [V(G + h) + v B(G + h)/c 2 Bt] , (1.37b) 

(G 23 , G 31 , G 12 ) = - TV (G + h) x v/c. (1.37c) 


Because ( G 23 , G 31 , G 12 ) v=0 = 0, we have for the tensors G jk and G jk that 
are analogous to the tensors F ik and F jk in the case of the electromagnetic 
field 


* and 


G jk = G jk 


(1.38a) 


G jk = 0. 


(1.38b) 
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II. LARMOR THEOREM FOR INVISCID FLOW 


In this section and the next one we shall neglect the effects of viscosity. 
Thus we assume that in (1.36) 

9 k SJ k □ 0 (2.1) 

where S jk is the viscosity stress tensor. 

As a preliminary to rewriting the fluid equation of motion (1.36), we note 
that the two terms of this equation that involve the variable particle mass y can 
be rewritten as follows: 

d (yv*)/dr - B^ (fic 2 ) = - [3J (yv k ) - B k (/i. W)] v k . (2.2) 

Making use of (1.5), (2.1), and (2.2), we find that (1.36) can be written in the fol- 
lowing form: 

[(B j p k - 3 k p j ) + (m/c) 0J k ] v k = 0 (2.3) 

where p' is the canonical particle momentum defined by the relation 


pJ = y v J + (q/c) A J = (p°, p); 

(2.4a) 

P° = (y* c 2 + q A°/c) , 

(2.4b) 

p = (y* v + q A/c) , 

(2.4c) 


where y* — Vy is the variable particle mass in the observer's frame. 
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We now introduce the antisymmetric canonical vorticity tensor V J k , defined 


as follows: 


V jk = - ( cb p k - d k pj ) ; 

(2.5a) 

(V 10 , V 20 , V 30 ) = Bp/cBt + Vp° = F/c 

(2.5b) 

(V 23 , V 31 , V 12 ) = V x p = V; 

(2.5c) 


where (2.5b) and (2.5c) constitute the definitions of the 3-force F and the 3- 
vorticity V in terms of the components of V ik . 

We note now that (2.3) implies that 

W k = (m/c) 0 jk + 2/U. ib k (2.6) 

where v jk is an antisymmetric tensor that is undetermined except for the 
condition 

^ jk v k = 0 (2.7) 


which implies that it can be written in terms of a completely arbitrary 3-vector 
v as follows: 


(v 10 , i/ 20 t y 30 ) = - Tv x v/c; 


( [v 23 , v 31 , v 12 ) = T v 


(2.8a) 

(2.8b) 


We shall see that v iV can be interpreted as the tensor that describes the residual 
rotation of the fluid, i.e. that part of the total rotation Q 1 k that is not produced 
by the action of the external fields on the fluid, but rather must be regarded as 
one of the initial starting conditions of the fluid which has been preserved 
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unchanged. (It was to enhance this interpretation that in (2.6) the factor 2 /a was 
inserted in front of ^ jk .) 

Using (1.32), (2.5), and (2.8), we can rewrite (2.6) as the following two 3- 
vector equations: 

F = 3p/3t + c V p° = m*(©_ 0xv/c)-2m*vxv; (2.9a) 

V = (m*/c) (0 + ©x v/c) + 2/U* v . (2.9b) 

The first of these is an alternative form of the equation of motion. The second 
expresses the canonical 3-vorticity V in terms of the residual rotation v and 
the effects of the entropy-dependent forces. 

Using (1.16a), (1.37a), and (2.4a) in the definition of V jk given in (2.5a), we 
arrive at the following relation: 

V jk = 2/u.Oi k _ ( m / c )G jk - (q/c)F jk . (2.10) 

Making use of this to eliminate V jk in (2.6), we arrive at the desired relativistic 
generalization of the Larmor theorem for inviscid flow of a charged fluid: 

U jk = (q/2 fj. c) F j k + (m/2 /xc) G j k + (m/2 /xc) k + v iV . (2.11) 

The first term on the right is the covariant representation of the familiar Larmor 
precession. The second and third terms are completely analogous precessions 
that are produced by the /^-dependent and entropy-dependent forces respectively. 
(By " //-dependent force" we mean, of course, the sum of the gravitational and 
enthalpy-dependent forces.) The last term on the right side is completely un- 
determined except for the condition (2.7). Obviously this represents the fluid 
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rotation that would be present in the absence of any forces, and so must be re- 
garded as an initial starting condition of the fluid. Note that the tensors G jk and 
0l k have been defined in such a way that their coefficients in (2.11) have the same 
form as that of F ik except that the charge q is replaced by the mass m. 

The physical significance of (2.11) becomes most apparent if it is decomposed 
into one equation in which all the space-time components vanish in the fluid rest- 
frame, and a second equation in which all the space-space components vanish in 
this frame. The first of these equations is 

U jk = (q/2Mc)F jk + (m/2Mc)0 jk + v jk , (2.12) 


and the second is 

n jk = (q/2/xc) F )k + (m/2/xc)G jk + (m/2/uc) 0 jk . (2.13) 

Using (1.10c), (1.21b), (1.31c), and (2.8b), we find that the space-space part 
of (2.12) is equivalent to the following 3-vector statement of the Larmor theorem: 

Q = - (q/2/U.c) B + (m/2/xc) 0 + v (2.14) 

where B is the effective magnetic field defined in (l.lOd) and 0 is the analogous 
thermal field defined in (1.31d). Note that the ^-dependent force, which is a 
potential force, produces no precession in the fluid. As indicated in (1.22), the 
total rotation Q differs from the Q given above only by the Thomas precession 
Q t . The 3-vector equation that corresponds to the space-time components of 
(2.12) is obtained simply by taking the vector cross-product of (2.14) with the 
fluid 3-velocity v . 
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Using (1.9b), (1.18b), (1.28c), and (1.37b), we find that the space-time part 
of (2.13) is equivalent to the following 3-vector equation: 

p* dv/d r = - m V (G + h) + qE + m0-vd [m(G + h)/c 2 ] /Bt , (2.15) 

where E is the effective electric field intensity defined in (1.9d) and 0 is the 
analogous entropy force defined in (1.28e). This is an alternative form of the 
fluid equation of motion that is characterized by the fact that the particle mass 
p*, which includes the speed-dependent mass increase, stands outside the time 
derivative on the left. 

We can cast this equation of motion into a more explicit form by making 
use of (1.9d), (1.28e), and (5.2) of I: 

p* p* Dv = - pmVG - VP + p*q [(E-Bxv/c) - (E • v) v/c 2 ] 

- v [pm9 G/B t + 3P/3t + (p*T*mDs/c 2 )w-v]/c 2 

+ (p’T’mDs/c 2 ) w. (2.16) 

The last term on the right side of (2.16) is a force that arises when heat 
energy is transferred from the reservoir to the fluid. This energy has a mass 
associated with it, and the time rate at which heat mass is injected into unit 
volume of the fluid is p* T*mDs/c 2 . Because this mass has the velocity w rela- 
tive to the fluid, its absorption produces a drag force (p* T *mDs/c 2 ) w that tends 
to accelerate the fluid in the direction of the reservoir velocity. 
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As was pointed out in Section IV of I, this drag force performs the work 
(p* T*mDs/c 2 ) w • v per unit volume of the fluid. Work is also performed on the 
fluid by increases in the gravitational potential G and the pressure P. (cf. (2.18)of I). 
These three forms of work, unlike work performed by the electromagnetic field, cause 
an increase in the rest-mass of the fluid. Such an increase in rest-mass, with- 
out the application of a corresponding force to accelerate the newly added mass 
to the fluid velocity, causes the fluid to decelerate so as to conserve its momen- 
tum. This effect results in a drag force in the direction -v. This is the explana- 
tion of the next- to- last term on the right side of (2.16). 

We have shown that the Larmor theorem (2.14) and the equation of motion 
(2.16) are parts of a single tensor equation given in (2.11) that may be regarded 
as the covariant statement of the Larmor theorem. The contraction of (2.11) 
with the fluid 4-velocity yields the 4-vector equation of motion given in (1.36) 

(for 3 k S jk = 0). We shall now apply this formalism to derive the relativistic 
generalization of the Helmholtz equation for inviscid flow of a charged fluid in 
the presence of electromagnetic and gravitational fields. 
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m. HELMHOLTZ THEOREM FOR INVISCID FLOW 


As a preliminary to deriving the generalized Helmholtz equation, we must 
first derive an alternative to the expression for F given in (2.9a). If we define 
the 3-vector o' as 

S' =S+6xy/c, (3.1) 

then it follows from (1.25c) and (1.32) that 

(0 23 , 0 31 , 0 12 ) = r§' = F[(T V s) x (v + w) /c]- (3.2a) 

With this we can derive the following alternative to (1.25b): 

(0 10 , 0 20 , 0 30 ) = - r§' X v/c + (T V s)/r + (Ts/c 2 ) (v + w). (3.2b) 

From (2.5), (2.6), (2.8), and (3.2) we find 

F = c(V 10 , V 20 , V 30 ) = v x V + (m T Vs)/T + (mT s/c 2 )(v + w) (3.3a) 

where 

V = (V 23 , V 31 , V 12 ) = (m* /c)0'+ 2 m* v (3.3b) 

Equation (3.3a) is the desired expression for F. 

The covariant vorticity conservation equation is 

d’ V k ^ + 3 k V^ j + B^V jk = 0, (3.4a) 

which follows immediately from the form of the definition of V jk given in (2.5a). 
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(3.4a) is equivalent to the following two equations: 


V • V = 0; (3.4b) 

BV/Bt = Vx F- (3.4c) 

Using (3.3a) in (3.4c) we have 

BV/Bt = Vx (v x V) + Vx [(mT V s )/ F ] + Vx [(mTs/c 2 )(v + w) ] . (3.5) 

Substituting (3.4b) and the continuity equation 

V • V = - (D p* ) / /O* (3.6) 

into the identity 

Vx(v X V) = - v • V V + v (V ■ V) t V • Vv - V( V • v) (3-7) 

and using the resulting relation in (3.5), we arrive at the generalized Helmholtz 

equation: 

D( S/p*) = {S/p*) ■ Vv + (l/p* ) V x [(i _ v 2 /c 2 ) m T* Vs ] 

+(l/p*) Vx [ (m T* Ds /c 2 ) (v + w) ] , (3.8) 

where in the last two terms we have made use of the fact that T./F = (1 - v 2 /c 2 )T* 
and Ts = T*Ds. 

Because 

V=Vxp = Vx( / u*v + q A/ c) 

Vx (u*v) + (q/c) B (3.9) 
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where 


p* =m(l - v 2 /c 2 )" 1/2 (1 + G/c 2 + h/c 2 ), (3.10) 

the effects of the electromagnetic, gravitational, and thermal fields are included 
in (3.8), as well as all relativistic effects. 

The last term on the right side of (3.8) arises from the injection of heat into 
the fluid and vanishes when Ds= 0. The next-to-last term vanishes if the specific 
entropy is everywhere the same. Thus if the flow is isentropic (Ds = Vs = 0), the 
relativistic Helmholtz equation for inviscid flow of a charged fluid in the presence 
of electromagnetic, gravitational, and thermal (enthalpy) fields reduces to the 
familiar form D (V / p*) = (V/ p*) • Vv if the generalized vorticity is defined as in 
(3. 9). 5 

(except in the ratio q/c) . 

The nonrelativistic limit is obtained by letting c go to infinity Then p*- p, 
T* T , and (3.8) becomes 

D[(Vxv + q B/mc)/,o] = [(V x v + qB/mc)/p] •Vv +(VT x Vs) /p. (3.11) 

(nonrelativistic) 

For B = 0 this is just the well-known Vazsonyi 6 generalization of the Helmholtz 
equation for the case of nonisentropic flow of an uncharged gas. 

If we go to the limit m - 0 in (3.8), we have 

(m = 0) D(B/p*)= (B/p*) • V v . (3.12) 
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This is the equation that expresses the fact that the magnetic field is "frozen 
into" a perfectly conducting fluid which, of course, is what we would have if the 
charge-carrying particles had zero mass. 

Having found a generalization of the Helmholtz equation that comes so close 
to the simple form D (V/ p*) = ( v/p*) • V v , we now pose the question whether it is 
possible to modify the definition of the vorticity V in such a way that this simple 
form of the equation will hold exactly . In the next section it will be shown that 
this is indeed possible, even when viscosity must be taken into account. The 
procedure consists of lumping together the entropy-dependent forces and the 
viscous forces and describing them in terms of a 4-vector function a j whose 
dynamical effects are completely analogous to those of the electromagnetic 4- 
potential A' . When the definition of the canonical particle momentum is modified 
to include a contribution from a 1 , we shall find that the vorticity, defined as the 
curl of this generalized canonical momentum, satisfies the simple form of Helm- 
holtz's equation exactly. 
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IV. GENERALIZED HELMHOLTZ THEOREM 


Using (5.13) of I in (1.36) above, the fluid equation of motion including vis- 
cosity ( B S j k ¥ 0) becomes 

d(/zv^)/dr = B* (mc 2 ) +(q/c)F^ k v k + tt* -mTB^s + 77 *, (4.1) 

where v > is the time rate of energy-momentum injection resulting from heat 
injection which, using the equations (2.20) and (5.13) of I, can be written 

77 J = -(B k Q jk )//? = (m/c) 0 jk v k + rnTB* s. (4.2) 

V 1 =-(\S*)/p (4.3) 

is the viscous force. The symmetric tensors Q jk and S ik are, of course, the 
heat and viscous stress contributions respectively to the fluid stress-energy 
tensor. 

Using (1.5) and (2.2), we can cast (4.1) into the following form: 

- [ B* (/xv k + qA k /c) - & (yuv^ + q A* / c) ] v fc = 7T> -m T B^ s + 77 J . (4.4) 

It is obvious that the contraction of the left side of this equation with v. vanishes. 
Thus we conclude that 


v . (jt ] - m T B^ s + 7 )> ) = 0. (4.5) 

Because the force (tt> - mT B j s+ r]> ) must satisfy this orthogonality condition, 
it can be expressed as the contraction of a suitably chosen antisymmetric tensor 
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with v k . Thus we define the antisymmetric tensor f jk by means of the follow- 
ing relation: 

(m/c)f jk v k - tt 1 - ml s -t r; j . (4.6a) 

Because 

f i k = _ f k i , (4.6b) 


it follows that the orthogonality condition (4.5) is automatically satisfied. Re- 
ferring to (4.2) and (4.3), it is evident that the defining equation for f jk can also 
be written in either of the following forms: 

fi k v k = 0> k v k + (c/m) v 1 ; (4.6c) 

or 

f jk v k = - c [TcP s + (1/pm) B k (Q jk + S jk )] . (4.6d) 

Because of the analogy between f ik and the electromagnetic field tensor F jk , 
the space-time components of f will be designated as e and the space-space 
components as b : 


(f 10 , 

f 20 , f 30 ) = e ; 

(4.7a) 

(f 23 . 

f 31 , f 12 ) = - b. 

(4.7b) 


The equation (4.6a) specifies only e , the value of e in the fluid rest frame. 

It says nothing about b . We can remove this indeterminancy by requiring that 
f ik be the curl of a 4-vector : 

fi k = a k - ^ a> . (4.8) 
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Inasmuch as the addition to ai of the 4-gradient of an arbitrary scalar function 
has no effect on f jk , the 4-vector a 1 has only three significant degrees of free- 
dom, which corresponds to the fact that the 4-force (tt> - mT s + p j ) has three 
degrees of freedom. (One degree of freedom is removed by the orthogonality 
requirement (4.5).) The equation of motion for a j follows from (4.6a) and (4.8): 

da j /dr = v k 3* a k - (c/m) (jt' s -idT3' s + rf ). (4.9) 

By introducing f jk and we have placed the handling of the viscous and 
entropy-dependent forces on the same footing as the handling of the electromag- 
netic forces. In the case of the latter, however, the appropriate field equations 
are Maxwells equations, whereas in the case of the former, the "field equations" 
are given by either (4.6) or (4.9). This point of view will be somewhat more 
fully developed in the following paper. 

Our present purpose, however, is to derive the generalized Helmholtz equa- 
tion. To this end we note that if we substitute (4.6a) into (4.4), the resulting 
equation of motion is just the contraction of v k with the following equation: 

V jk + (m/c)fi k = 2/W k , (4.10) 

where V jk is the canonical vorticity defined in (2.5), and c <d k is an arbitrary 
antisymmetric tensor that satisfies the orthogonality condition 

ad k v k = 0 (4.11) 

which implies that ad k has the form 

(ai 10 , c<j 20 , oj 30 ) = — r a> x v /c, (4.12a) 
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(a) 23 , OJ 31 , 


w 12 ) = Tcy 


(4.12b) 


where (4.12b) is the definition of the 3-vector in terms of the components of 

The tensor 2p.o)> k will play the role of a generalized vorticity, which we shall 
call the intrinsic vorticity*. It can be written in the form 

2 /jLC0> k = (3* j» k - (4.13) 

where 

= n v* + (q/c)A^ + (m/c) a’ (4.14) 

is the generalized canonical particle momentum that includes the contribution 
(m/c) a’ that arises from the viscous and entropy-dependent forces. 

Substituting (2.10) into (4.10), we arrive at the generalized Larmor theorem: 
fii k = (q/2Mc)F jk + (m/2/xc)G jk + (m/2/ic)f ik + cJ k . (4.15) 

A comparison of this with (2.11) shows that 

(m/ 2Mc)(0 jk - f jk ) = ct> ik - i> jk . (4.16) 

Contracting this with v k and using (2.7) and (4.11), we have 

(0i k _ fi k )v k = 0. (4.17) 

Thus f i k differs from 0i k at most by a tensor that is orthogonal to v fc . 

The tensor v' k was interpreted as the residual rotation, i.e. the rotation 
not produced by external forces, but rather associated with the initial conditions 
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of the fluid. From (4.15) it is evident that a similar interpretation may be as- 


signed to od k . To distinguish it from v> k , we shall call a> >k the intrinsic (rather 
than the residual ) rotation. Of these two tensors, which in general are not equal, 
o4 k is the more fundamental. This follows from the fact that it is the flux of &> 

(or more exactly, the flux of the intrinsic vorticity 2 and not that of v that 
in general is conserved. We shall now prove this. 

From (4.13) it follows directly that 

3>(Ai<u k/t ) + 3 k ( M «/i) + 3- (, iiwi k ) = 0, (4.18a) 

which is equivalent to the following two equations: 

V-( M *o>) = 0, (4.18b) 

d(/u* &>)/3t = Vx (V X ju* <y ), (4.18c) 

where use has been made of (4.12). It is well-known 7 that these two equations 
imply the validity of the following conservation law: 

D J n*o) • dS = 0 (4.19) 

s 

where D = (d/dt) + v • V is the substantial time differentiation operator, and the 
integral is taken over any surface S which moves -with the fluid. Thus the flux 
of the intrinsic vorticity 2 p*<w passing through any closed loop that moves with 
the fluid never diffuses out of this enclosed area. It is easily verified that this 
would be true for the 2/jl*v flux only if 0 jk , like f jk , were expressible as the 
curl of a 4-vector. 
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Because (4.18c) lacks the two entropy-dependent terms that appear in (3.5), 
when we duplicate the steps that led from (3.5) to (3.8) starting now from (4.18e), 
we arrive at the following simple form for the generalized Helmholtz equation: 

D (p*«/p*) - 0**>/p*) • Vv. (4.20) 

(The ratio p*/ p* could of course also be written p/p, but it seems more appro- 
priate to refer everything to the observer’s frame.) 

It is possible to give this generalized Helmholtz equation an interesting 
microscopic interpretation which will now be discussed. 

Microscopic Interpretation of Helmholtz Equation 

We may take the point of view that the mass p* of each of the p* particles 
that occupy unit volume of the fluid is uniformly distributed over a volume Uj 
= 1 / (J* which for the sake of simplicity we may consider to have the form of a 
right cylinder of height -£ and cross-sectional area Q. We assume that at some 
given instant of time the axes of the particle-cylinders are all parallel to the 
local value of the intrinsic rotation o>, and that the cylinders are rotating about 
their axes with the angular velocity «. The total angular momentum of any 
particle may be regarded as the sum of its external and internal parts. The 
external part is just the angular momentum about some specified origin of a 
point mass p* located at the center of the particle-cylinder. It thus depends 
only on the position and displacement velocity of the cylinder, but not on its 
rotation about its own axis. The internal angular momentum a is just the product 
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of the angular velocity « of the cylinder and its moment of inertia which is 

*Cx / 2tt = p* /(2vlp*): 

a = /A>/(2 tt ip*). (4.21) 

We further stipulate that the ends of each particle-cylinder be fixed in the 
fluid. Thus if £ is the vector connecting the two ends of the cylinder, it must 
obey the well-known equation 8 for the time dependence of the vector displace- 
ment between two neighboring fluid particles: 

D£ = £ -V v . (4.22) 

But this equation has the same form as (4.20). This implies that 

M*»/ p* = K£ (4.23) 

where K is a constant along the trajectory of any given particle, i.e. DK = 0. 

Thus if at any instant of time the axis of a particle-cylinder is parallel to the 
local value of w, then as the cylinder is swept along with the fluid it will find 
that its orientation is always such that its axis is parallel to the local value of 
cj. Although the direction of a changes in general as the cylinder is carried along 
with the fluid, its magnitude |<r| is constant. This follows from (4.23): 

a = p.*o>/(2n-ip*) = (K '277-) (4.24) 

where £ is the unit vector along the cylinder axis. Because DK = 0, we have 
then 


D|o| = D(K/2 tt) = 0. 


(4.25) 
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The behavior of a particle-cylinder whose two ends are fixed to the fluid 
points A and B is illustrated in Figure 1. As the particle-cylinder is carried 
along with the fluid, it changes (in general) its size and shape, and hence its 
moment of inertia, but the magnitude of its internal angular momentum a = joj 
remains constant, which means that the moment of inertia must vary inversely 
with [ co ! • 

Note that the internal angular velocity of the particle-cylinders is the in- 
trinsic angular velocity co and not the fluid rotation Q . These two differ by the 
precession produced by the external fields. We conclude then that the external 
fields are unable to apply a torque about the axis of any given particle-cylinder. 
Because of this its angular momentum a is conserved except for the fact that, 
because the cylinder is embedded in the fluid, the axis of the cylinder, and hence 
the direction of o, changes with time. Thus, if the particle-cylinders do not have 
any internal angular momentum, there is no way in which they can acquire it. 
This, of course, corresponds to the well-known implication of the Helmholtz 
equation that, if the intrinsic vorticity 2p*« vanishes at any point on a particle 
trajectory, it must be zero along the entire trajectory. 

It would be of interest to derive the equation of motion of a . From (4.20) 
and (4.21) 


Do = a -Vv - o(D l)/l. 


(4.26) 


Using (4.22) we find 


(D4)A£ = C 2 (Df ) • £ =(£-V v ).£/^2 


- (o '^7 v) • a /cr 2 . 


(4.27) 
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Thus (4.26) becomes 


Da = a ■ Vv - a(a-Vv)-a/cr 2 , 


(4.28a) 


or 


Da = ax [ (a • V v ) x a ] ,/cr 2 . 


(4.28b) 


The second of these two equations shows that the motion of a may be regarded 
as a precession with the a -dependent angular velocity -(a • V v ) xa/o- 2 . 

The dependence of the motion of a on the fluid velocity can be made more 
explicit by rewriting (4.28a) in the following form: 


Da 


[2) - a - a(a , ®'a)/cr 2 ] ■+ — (Vxv)xa, 


where 2) is the symmetric fluid distortion dyadic given by 


(4.29a) 


© =— [Vv + (V v ) ]. 

2 y c 


(4.29b) 


where (Vv ) is the conjugate (or transpose) of Vv. The second term on the right 
side of (4.29a) is simply a precession of angular velocity — (V x v) , which is 
the nonrelativistic fluid rotation. Thus this term is simply the precession of the 
internal angular momentum that results because the axis of the particle-cylinder is 
imbedded in the fluid and must rotate with it. The first term on the right side of 
(4.29a) describes the change in the orientation of the cylinder axis (and hence 
of a ) that results from fluid distortion. 
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The precession — V x v may also be expressed in terms of the external 
fields. Using (1.4b), (1.16c), (1.37c), (4.7b), and (4.12b) in (4.15) we find 


-L V x v = - (2p*c) 1 {qB + [mV (G + h) + /xVF c 2 ] x v /c + mb } + w. (4.30) 


Because « and a are parallel, « x a = 0. Thus the precession velocity of a in 
(4.29a) is given by the first term on the right side of (4.30). The middle term 
inside the braces is a relativistic precession which results from the gravita- 
tional, thermal, and Bernoulli forces, and which vanishes in the rest frame of 
the fluid. 

Equation (4.21) can be used to estimate the magnitude a of the internal 
angular momentum, if we assume that l ^ Dj /3 = ( p*) ~ 1/3 . Taking i* to be of 
the order of the proton mass, we find that for laboratory and astrophysical 
problems 


0 < cr <, 10~ 36 er g- sec . (4.31) 

It is interesting to note that if we were to take a classical hydrodynamic al view 
of a large nucleus, assuming that o> ~ 0.1 c/R where the nuclear radius R is of 
the order of 10* 13 cm, then we would find that cr ~ 10" 27 erg-sec, where +i is 
Planck's constant. The angular momenta that one encounters in nuclear and 
atomic problems are, of course, of the order of magnitude of "ft. Thus the in- 
ternal angular momentum a that we encounter in classical problems is much 
smaller than the typical magnitude of the angular momentum per particle that 
we encounter in atomic and nuclear problems, and would approach this magni- 
tude only if we were to attempt to apply the classical formalism to problems 
of atomic and nuclear dimensions. 
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V. CONCLUSIONS 


The formalism based on the description of the viscous and entropy-dependent 
forces in terms of a 4-potential a j that is analogous in its dynamical effects to 
the electromagnetic 4-potential has two advantages over the formalism based on 
the entropy force tensor 0 jk : First, the 4-potential a 1 can be used to describe 
the viscous as well as the entropy-dependent forces. Second, the intrinsic vor- 
ticity 2 J u*co in the a’ -formalism is more fundamental than the canonical vorticity 
V x (ju*v + qA /c) that enters into the 0 ik -formalism in that the Helmholtz equa- 
tion in the former case yields an exact conservation law for the 2p*co flux 
whereas the canonical vorticity flux in the latter case is not in general conserved. 
The most interesting feature of the a’ -formalism, however, is that it allows the 
fluid equation of motion to be cast into the form of a generalized Hamilton-Jacobi 
equation, whose derivation is given in the following paper. 
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